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Conventional approach to quantum mechanics in phase space, {q,p), is to take the operator based 
quantum mechanics of Schrodinger, or and equivalent, and assign a c-number function in phase 
space to it. We propose to begin with a higher level of abstraction, in which the independence and 
the symmetric role of q and p is maintained throughout, and at once arrive at phase space state 
functions. Upon reduction to the q- or p-space the proposed formalism gives the conventional quan- 
tum mechanics, however, with a definite rule for ordering of factors of non commuting observables. 
Further conceptual and practical merits of the formalism are demonstrated throughout the text. 

PACS numbers: 



I. INTRODUCTION 



Wigner's 1932 initiative is a reformulation of the operator based quantum theory of Schrodinger in the language 
of c-number distribution functions in a phase space. His prescription, however, turns out to have a feature extra to 
what one finds in Schrodinger 's theory. There is nothing in the founding principles of the operator based theory to 
prescribe a rule for ordering of the factors of non-commuting operators in a product. In contrast, Wigner's formalism, 
upon reduction from phase space to the configuration space, acquires Weyl's ordering 1^,1^. How and at what stage, 
in going from Schrodinger 's state functions in configuration space to those of Wigner in phase space and again coming 
back to the configuration space, Weyl's ordering creeps in? This feature is not imique to Wigner's functions. Other 
distributions exist in the literature, e.g., Kirkwood Husimi Margenau and Hill etc. Each of them carries 
its own ordering rule, with no precedence in the configuration space formalism. Can one conjecture that the phase 
space formulations of quantum mechanics are more complete than their configuration space counterpart, because of 
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their built-in ordering rules? If so, there should be a way to arrive at phase space formulations without reference 
to the conventional operator based theory. Here we argue that in the classical dynamics and classical statistical 
dynamics (Liouville's equation) the generalized coordinates and momenta, q and p, respectively, play symmetric and 
more importantly, independent roles. In the operator based quantum theory one or the other looses its identity at the 
expense of the other and the formalism reduces, to one in either q or p space. One could avoid this by carrying the q 
and p formalisms concomitantly and at once arrive at state functions in qp spaces. The so-obtained state functions 
are the qp representation of the mixed states of quantum statistical mechanics. The operator based theory emerges 
as a special case of this general one, but this time with a definite ordering rule for non-commutative operators. The 
rule depends on the nature of the q and p variables, adopted initially. 



II. EXTENSION OF THE CLASSICAL DYNAMICS 



Let q — {qi{t),i — 1,..,N} be the collection of the generalized coordinates describing the state of motion of a 
dynamical system. It is customary to assign a lagrangian, L'^{q,q), to the system, define the conjugate momenta, 
p = dL"^ /dq, and construct the H{q,p) = qp — L'^. One may do the other way around. Begin with a given H{q,p) and 
find L'^{q, q) as a solution of the following differential equation, 

One may, however, carry out the same procedure with q replaced hy p and arrive at a U'{p^p) satisfying the differential 
equation 

The use of to study the evolution of a dynamical system is not a common practice. But it is a possibility and has 
precedence Q. There is no barring to employ the two alternatives simultaneously. We follow Sobouti and Nasiri 
(hereafter, paper I) and define the "extended lagrangian" 

i^{q,q;p,p) = -qp-qp + L'^iq^q) + lp{p,p). (3) 

The first two terms on the right hand side constitute a total time derivative and are introduced for later convenience. 
One may now write down the Euler-Lagrange equations for q and p, 

dt dq dq dt dq dq 
d_dC _ dC _ d_dIJ^ _ dU^ _ ^ 
dt dp dp dt dp dp 

Equation l|4a|l is the conventional equation of motion in q space. With preassigned initial values <z(to) and q{tQ) at 
to it can be solved for the orbits q{t) in q space. Similarly, with given initial values p{to)a,nd p{to), Eq (|4b|l can be 
solved for the orbits p{t) in the p space. The condition for q and p orbits to represent the same state of motion of 



the system are p(<o) = 9L'' /dq\to and (/(to) = dL^/dpltg. Such a state of motion will be referred to as a "pure state". 
Otherwise it will be called a "mixed state" of motion. The nomenclature is from the statistical quantum mechanics 
and it will be seen later that they imply the same notions as therein. On a pure state p and q are initially canonically 
conjugate pairs and it is shown in paper I that once they are canonically conjugate at one time they remain so for all 
times. On the other hand there are no restrictions on the initial values of q and p on mixed states. Therefore, q{t) 
and p{t) remain unrelated and evolve independently. The existence of the extended lagrangian C{q, q;p,p)^ however, 
permits the following "extended momenta" to be defined, 

dC dLt 



dC _ dLP 
dp dp 



These in turn allow an "extended hamiltonian" to be defined through the following Legendre transformation, 

H{q,TTq;p,7Tp) = qiTq + pTTp - C{q,q;p,p). (6) 

To eliminate q and p from H one substitutes a) for £ from Eq Q , b) for L'' and from Eqs ^ and Q and c) for 
dL'^/dq and dL^ /dp from Eqs H5a|l and (|5b|) . One arrives at 



H{q,TTg;p,Trp) = H{q,p + 'Kq) - H{q + Trp ,p) 



y- 

^ Ti l 



n=0 



d'^H „ d'^H 
7r„ 



(7) 



where the derivatives are to be evaluated at {q,p). We leave it to the reader to familiarize him/herself with TL 
by writing down four Hamilton's equations for q, -Kq, p and tt'^ . Here, the condition for pure state motions is 
^9(^0) — '^p{to) = 0, and once they are initially zero they remain so for all times. Then by Eqs (|5a|l and Ij5b|l q and p 
turn into canonically conjugate pairs for all times (paper I). To summarize, for any dynamical system we introduce an 
extended phase space, (g, tt^; p, iTp), extended momenta, lagrangians, and hamiltonians. All concepts and procedures 
of the conventional dynamics are extendible to this extended dynamics. Of particular relevance to this paper, that 
will be referred to shortly, are: 1) canonical transformations from one set of variables (q, tt^; p, -Kp) to another, and 
2) Poisson's brackets extended as 

{F G} -——-—— + ——- —— (8) 
dq d-Kq duq dq dp d-Kp d-Kp dp 



III. QUANTUM DYNAMICS IN qp SPACE 

Now that we have the extended hamiltonian of Eq Q we may construct a quantum mechanics in qp space. We do 
this on the following premises: 



1. Let X be the function space of all integrable complex functions xiliP)- L^t q, TTg, p, and, tt^ be operators on 
X, satisfying the commutation rules 

[q, T^q] = [P, TTp] = ih, [q,p] = [tT,, TTp] = [q, TTp] = [p, TTg] = 0. (9) 

These are the fundamental Poisson brackets of Eq ijSJ, promoted to commutation brackets by Dirac's prescrip- 
tion. Note the manifest independence of q and p in the vanishing of their commutation brackets. 

2. By the virtue of Eq lO, Ti. is now an operator on X. Let xiliPjt) € X he a state function satisfying the 
Schrodinger-like equation 



d d 
H{q,p - ih—) - H{q - ^f^g^^P) 



X- (10) 

3. Let the rule to evaluate the expectation values of an observable 0{q,p), a real c-number operator on X, be 

< 0{q,p) J 0{q,p) Re xdqdp =^J 0(g,p)(x + X*)dqdp. (11) 

We will return to this averaging rule shortly, and revise it. The logic behind it, however, is to be noted, the 
averages of observables should be real. In what follows we demonstrate that, 1) the formalism so designed 
is a theory of quantum ensembles in phase space. Its pure state case is the conventional quantum mechanics, 
however, with a definite ordering rule accompanying it. 2) It can be transformed to other phase space formalisms, 
including to that of Wigner, by suitable unitary or similarity transformations on X . The latter in turn originates 
from suitable canonical transformations from one extended phase space coordinates to another. 



A. Solutions of Equation (IIUB 



To begin with, x is of the form 



x(g,p) = F(q,p)e"^P«/^ 



(12) 



The exponential factor is a consequence of the total time derivative, —d{qp)/dt in Eq It is easily verified that 

(13a) 



{p-ih—)x = ^n—e-^p'^/\ 



(q - ih^)x = 



dp 



OF 
dp 



-ipq/ h 



Substitution of Eqs (13) and (221 in Eq gives 



H[q^^^n^)-H{-^h^^,p) 



F. 



(13b) 



(14) 



The operators on the right hand side of Eq H14|l are recognized as the hamiltonians of the conventional quantum 
mechanics, the first in q and the second in p representation. Thus, one obtains the superposition of the separable 
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solutions 



where 



xiq,P,t)^J2^-^Ml^*WP^t)e-''"^\ (15) 

a, 13 



ih^^ = H{q,-ih -^)4>a, (16a) 



ih^=H{ih^, p^^. (16b) 



To each tpaio) there corresponds a 4>a{p) that are Fourier transforms of each other, 



^o.{q) = ttA^ I <f>a{p}e'^"'^^dp , (17) 



where N is the number of degrees of freedom of the system. 



B. The averaging rule revisited: Acceptable state functions 

Let Q{q) be an observable represented by a real polynomial or series in q. Its matrix representation, Q, in either 
X", ip-, or (/<-basis is hermitian. Thus 



Qpa = J Mq)Q{qWp{p)e-''"'^^dp dq 
= J i*p{q)Q{q)'>lJa{q)dq 

(l)l{p)Q{ih—)(j)^{p)dp = Qlp, (18) 

where we have used the fact that ip ^^nd (j) bases are the Fourier transforms of each other. The coefficient {2nh)^^^'^ 
is suppressed for brevity. The expectation value of Q, by Eq l|ll|l . now becomes 

<Q>^\j Q{X + Xldpdq = itr[Q(A + At)], (19) 

where A is the matrix of Aap of Eq (|15(l . This gives the freedom of choosing A = A"^ and of simplifying Eq Hll|l 
to read < Q >= f Qxdp dq =tr{QA). Choosing Q{q) = 1, imposes the further restriction trA = 1. Requiring the 
averages of all positive definite functions of q to be positive still restricts A to be a positive definite matrix. Had one 
chosen a differentiable function P{p) instead of Q{q), one still would have arrived at the same requirements for A. To 
summarize, x of Eq H15|) is a physically acceptable solution if 

A = A^ , positive definite, and trA — 1. (20) 

With this provision the averaging rule of Eq Hll|l for Q{q) + P{p) reduces to 

< Q{q) + P{p) j{Q + P)xdp dq. (21) 
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For a product Q{q)P{p), by the prescription of Eq and with the restrictions of Eq (|20|l on A, one has 



< QP Re tr(QPA) = - tr(QPA + APQ) = tr[-(QP + PQ)A], 



(22) 



where Q and P are the matrix representations of Q{q) and P{p) as in Eq H18|l . Translation of this to the q space 
language, say, is 



Thus, upon reduction of the formalism of the present paper to that of the q-space, the ordering rule associated with 
a product Q{q)P{p) is the symmetric ordering. It has emerged from the formalism itself, unlike the ad hoc ordering 
rules of the conventional quantum mechanics. 



It was stated earlier that the proposed dynamics is essentially that of the ensembles. Here we elaborate on this, and 
show that l)thc classical limit of the theory is Liouville's equation that governs the dynamics of classical ensembles. 
2) Its pure state case is Schrodinger's operator based theory. 3) In its full generality the theory gives von Neumann's 
density matrix and the evolution equation associated with it. 



In Eq (|1(J|I expanding the hamiltonian operators about {q,p), and retaining only the first terms in the expansion, 
gives 



This is the Liouville equation for the distribution function of classical ensembles. Its most general solutions are 
x[g(t),p(t)], where q{t) and p{t) are the classical trajectories in q and p spaces. The two trajectories may represent 
the same state of motion if they satisfy the conditions of initial canonical conjugacy narrated below Eqs (4) . Otherwise, 
they remain independent and evolve independently. It is this classical notion of independence that we have carried 
through to the quantum formalism. Let us also note that the reduction of the phase space evolution equation to the 
classical Liouville's equation is a common feature of all such formalisms. 




(23) 



IV. MORE ABOUT EQUATION ([Ttni 



A. Classical correspondence 




(24) 
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B. Schrodinger's case 

Allowance for only one term in Eq (|15|) reproduces the conventional quantum mechanics in all its details. Thus 



X = ^(g)0*(p)e-'P«/^ (25a) 
dtlj d 

ih^ = H{q,-ih—)^, (25b) 

ip and (j) Foureir transforms of each other, (25c) 

/, = = = (25d) 

< QWPlr) >=\j r\P('iliptQ{q) + Q{q)P{-ih^)]vdq. (25e) 



Heisenberg's uncertainty principle follows immediately from Eqs (25) that one may find in standard texts in quantum 
mechanics. The ordering rule of Eq H25e|l is, however, the added feature of the theory. 

C. Density matrix and Von Neumann's equation 

The state function of Eq H15|) , as it stands represents the state of an ensemble in a mixed state. If the matrix A 

is diagonalizcd to Aaj3 ~ AaSafi^ X reduces to x = ^ AaV'a^Qe"'^'^''- Upon integration over q or p one immediately 

recognizes A^ as the probability of the system to be in the state ^Jaili t) or 0q(p, t). One may however do better. Let 

{'0„(g)} be a complete orthonormal time independent basis set, and {0„(p)} be its Fourier replica. These basis sets 

are not required to be the eigenstates of H{q,p), though this is a possibility. Hereafter, to avoid the ambiguity, we use 

the latin subscripts to denote the members of the basis set and reserve greek subscripts to denote the solutions of Eqs 

(|16a|l and (|16b|) . Expansion of x in these bases assumes the form xil^Pji) ~ ^mn{t)tpn{q)4'mip)^~^^'^^'^ ■ Substituting 

this form in Eq IjlUfl . multiplying the resulting equation by V^* ((/)(/)„ (p)e'^'?/'', and integrating over q and p gives 

dA .... 
ih— = [A,H], A^ A'' positive definite and trA = 1, (26) 

where A is the matrix of the expansion coefficients and H that of the H{q,p) in either x-, i'- or 0- basis. Equation l|26|l 
is von Neumann's equation for the evolution of the density matrix. As is known the case tT:{A^)—iT:A — 1 represents 
an ensemble in a pure state. If tr{A'^) < 1, the ensemble is in a mixed state. 

V. CANONICAL TRANSFORMATIONS 

All machinery of the canonical transformations from one extended coordinate system to another and their associated 
unitary or similarity transformations in the function space are available for a forage of deliberations. Except for a 
passing remark on the prospects of fuller uses of this approach at the end of this section, here we confine ourselves to 
one one-parameter family of transformations of which Wigner's state function emerges as a special case. Husimi's all 
positive distribution functions are also briefly mentioned. 
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Consider the infinitesimal transformations 

q^Q - SaUp, TTq = IIq; p = P - SaUq, TTp = lip, (27) 
The generator of the transformation is G = ■Kp'Kq. To this (and for a finite a) there corresponds the unitary operator 

[/„ = e^ ^e*''"^, C/tf/, = 1, (28) 

in the function space. Operating by Ua on a pure state function xiq^Pi t) — ''P{q)4'*{p) Gxp{—ipq/h) generates another 
state function (let us call it a-representation) 

Xo(Q,p,t) = f/aX=(^)'^ / i^iq-aT)riq + {l-cx)r)e'P-/^dT. (29) 



See Appendix for proof of Eq For a = 1/2, Eq (|^ gives Wigner's standard function IsHlfll, Xi/2 = W{q,p,t). 
The cases a = and 1 simply give back x smd x* of this paper, respectively. Similarly, operation by Ua on Eq I|1U|) 
gives the evolution equation for Xa 

^^ dx.{q,P,t) ^ ^^d_^^UaX) = {U.nuDUaX, 

.^ dxc.{q,P,t) ff{l-2a) pd 9" 

dt ^ ^"^"^ ^^d^2X.-^f^-g-^X. + T. -I H^) ^d^X.. (30) 

See Appendix for proof of Eq 13UI) . For a ~ 1/2, even n terms in Eq (|30|) cancel out and one again recovers Wigner's 
evolution equation, (^j. See Eq (A7). 



A. Assigning g-space operators to phase space functions; ordering rule 

The phase space state functions are devised to evaluate the expectation values of a c- number observable, F(q,p), 
by integrations over the phase space. Upon reduction to the q space, say, /(q^p) turns into a differential operator 
in terms of q and tt^. The questions are: l)how different factors of non commuting q and tt^ are ordered in a given 
a-representation? 2) Averaging a given F{q,p) with different Xa's gives different values, how such averages change 
from one a-representation to another? Let Fa{q, i^q) be the q space operator corresponding to the c- number monomial 
qTipm phase space when averaged by Xa- The defining equation for Fa{q, tt^) is 

< j qy'^Xcdpdq = J r{q)Fa{q,7rg)tP{q)dq. (31) 

For the combination of a = and 1 corresponding to {x + X*) of Eq Hll|) this is already worked out in Eq H25e|l and 
is the symmetric ordering 

«>"^^(«"< (32) 



For a general a, it is given in Appendix Eq (AlO) 



g>™ ^ E ((1 - «)T,)'-9"(a7r,)™-^ (33) 



m 

r=o \ r 



For a — 1/2 this reduces to Weyl's ordering |a| which is known to go with Wigner's functions. To answer the 
second question we note the foUowing 

< g>" >a=o^ I q^p'^xdq dp^ I q''p"'Ulxc.dq dp ^ J U^iq^^Xadq dp =< C/„((7>") >„, (34) 

where by Eq (|29|l we have used x — U^Xa- The conclusion is that g^p™ averaged by x is the same as Ua{q"'p™) 
averaged by Xa ■ Upon adoption of Ua ^ J2 ^~^k\^ ^ '§pk and operation by it on one finds 



TO 



smaller of n or ni 



-k 



qnpm ^ (_iftQ,)„^nq"~lp™-l + l{-ihafn{n - f)TO(m - l)q"-^p"'~^ + ... (35) 



B. Assigning phase space functions to q space operators 



To a given operator F{q,T:q), a Taylor-expanded series in whatever order of powers of q and TTg, we associate the 
following c-nuniber function 



F{q,P) = X! F^nXnrn q\F\p > e~ 



(36) 



where F,„„ =< n\F\m >, is the matrix element of F in the basis of the eigenstates of H{q, vr^). This is actually the 
inverse of the procedure that we used in Eq H31|) to associate an operator with a c-number function (let a = and 
replace g"p™ by F{q,p) in Eq H33|l to see the analogy). The second equality in Eq (|35|l expresses the same in the ket- 
and bra- notation of Dirac. 

The corresponding function in a-representation is simply 



F^{q,p) = U^F{q,p) = E^^™„!7„Xnm = / < « " "^l^^ + " > e'^^'^^dr. 



(37) 



This is actually the generalization of Eq H29|l for a general operator F{q,nq). The rule for the product F = AB is 
worked out in Eq (A13): 



n=0 



Fiq,p=) < q\A{q,TTq)B{q,7Tq)\p > e'^^'^l^ = ^ 
One may also work out the a- representation of Eq (|38|l : 
F„(g,p) = UaF{q,p) = A„ 



dp"" 



dq" 



(38) 



q + tna—,p — ina(l — a) — 
op oq 



E 

n—O 



n'. 



_d d_ 

dqA dpB 



{1-a) 



Baiq,p) 
d d 



dp_ 



A oqs 



Aa{q,p)Ba {q,p). 



(39) 
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See Eq (A17) for details of the derivation. In section VI, we analyze Bloch's problem as an illustration of the use of 
the developments of the last two subsections. 



C. A remark on general transformations: 

An economical way of treating canonical transformations is the symplectic formalism. Let ij be the column vector 
(q, TTq] p, TTp). The equations of the classical dynamics assume the following form 

^-jf, (40) 
where 'H{ri) is the extended hamiltonian of Eq ITJ and J is the symplectic metric 

l\ 

(41) 

-1 oj 

An infinitesimal canonical transformation from t/ t n r/ I Sr/ is of the form 

fl + Sfl = V-e^^ , (42) 

where G is the generator of the transformation and e indicates its infinitesimal character. The matrix of the 
transformation is 

The condition for canonicity is 




= % - . (43) 



MJM^ = J + 0(e2). (44) 

This imposes the condition on G to be either linear in rjj or quadratic and symmetric in rji, rjj or both. For clarity, 
here after we confine our discussion to a system of one degree of freedom, = 1 . The most general form of G with 
the restriction just mentioned is 

G{ri) = Qi-qi + aij-qi-qj, i,j = 1, 2,3,4 correspond to q, TTq; p, TTp, (45) 

where the four parameters a; initiate translations and the 10 symmetric aij cause rotations, boosts, squeezes, scale 
changes etc. The ten transformations aij constitute a symplectic group SP{4), and is locally isomorphic to the (3 + 2)- 
dimensional Lorentz group. This is the group that Kim and Noz ll] encounter in their study of four dimensional 
phase space consisting of two oscillators representation of 0(3 + 2), and proves to be a useful! mathematical tool 
in quantum optics. To each of the fourteen transformations of Eq l|45|) there corresponds a unitary or similarity 
transformation in the function space. That of Eq (|27|) it is unitary. An example of non unitary operators is the 
following. To the canonical coordinate transformation 

ie 1 

q = Q + + -Hp, = Uq 



P^P+ ^Hp + ^Hq, TTp = Hp, (46) 
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there corresponds the complex similarity operator 

Se = exp 



£5^ f^d^. 
M9g2 + 4s dp^' 



ih 52 



2 dqdp 

where e is a finite parameter of the transformation. Husimi's 4] all positive distribution in terms of x is 



(47) 



XHus{q,P,e) = Sex{q,p)- 



(48) 



VI. BLOCH'S EQUATION IN PHASE SPACE 

In this section we intend to illustrate some usage of the formalism developed so far. In any discussion of statistical 
mechanics, the partition function, Z{(3)~trQ, ri(g, tt^) = exp(— /3i7), plays a pivotal role. Its calculation , ho wever, 
is often cumbersome. O ne p ractice is to translate Ct and the corresponding Bloch's differential equation '1^ into a 
phase space language [sllioll. solve the equation for a c-number il{q,p,/3) and calculate Z{f3) — J il{q,p, /3)dqdp. The 
ease of doing the job depends on the choice of the c-number assigned to ft. Our suggestion is that of Eq H3t)|l . Bloch's 
equation for fl is 

(90 

— = -Hn = -nH. (49) 

We apply the rule of Eq (|36|l to Eq 149|) . Noting that H{q^Tiq) H{q,p) — p^ /2m + V{q) and using the product 
rule of Eq lp!H|l gives 



I I dp [ ' m dq 2m dq^ 

This same result is obtained in 13], however, by a totally different approach and through much lengthier calculations 
using Moyal's characteristic technique. Equation H49|l in Wigner's representation is obtained by replacing Cl with 
^w{QiP] P), H with I2m + V{q) and using Eq (A17) with a = 1/2 to find the expression corresponding to HCi. In 
agreement with Q| one finds 

dp — = 1 2;^ - ^^d^^^'^w^] (51) 

The contrast between the two Eqs (|50|l and H51(l is striking. The former is a second order differential equation in 
q and the exact quantum effects in it appear as h and only, while the latter in addition to is an nth order 
differential equations in p and has all powers of /i in it. In the following we solve Eq H50() and give the partition 
functions for the simple harmonic and linear potentials. 



A. Simple harmonic potential 

For H{q,p) — l/2{p'^/m + m.ui'^q^) the solution is of the form 



= exp 



-A{p)H{q,p)-iB{p)^-C{p) . (52) 



Substituting this in Eq (|50|l and letting the coefRcients of different powers of q and p vanish, gives 

dA 



1 - h'uj'A^ 

h^u'^A{l - B), 
1 . 



dp 
dB 

^ = -^-hWA. 
dp 2 
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(53a) 
(53b) 
(53c) 



The condition n{q,p,0) = 1 imposes the boundary conditions A{0) — B{0) — C(0) = 0. With these provisions one 
finds 



A{P) = -^tanh Pfuv, 
nw 

B{P) = tanh /3?kj tanh 
C{P) = -ilncosh/3fta;. 



Phoj 



The partition function is 



Z{P)^ 



2Trh 



-\-oo 



n{q,p,P)dqdp 



2 sinh 



phuj 



1 -1 



(54a) 
(54b) 
(54c) 

(55) 



The normalized density function is xiQiP^/^) — ^(^iP? /3)/27rfiZ(/?), with low and high temperature limits 

X — exp[-H{q,p) - — ] phw > 1 
nn a 



Plu 
2^ 



exp [-PH{q,p)] 



Phw < 1, 



(56) 



in agreement with the quantum and classical limits, respectively. 



B. Linear potential 

The case is of interest for quark model where a sea of semi infinite matter creates a linear potential V{q) = kq, 
< q < oo, and fc > 0. By the same procedure above one obtains 

ip^phk p^h^k^' 



^{q,P,P) = exp 



-PH 



2m 



6m 



(57) 



ZiP) 



27rm 



exp 



P^feP_ 
24m 



(58) 



27rm 







exp 



-PH 



iP'^-phk P^fi^k^ 



2m 



8m 



The corresponding Wigner's function can be obtained by letting U1/2 operate on Eq H59|l . 



(59) 
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VII. CONCLUSION 



We have developed a quantum mechanics in phase space by carrying the independent and symmetric roles of q 
and p, so eminent in the hamiltonian formulation of the classical mechanics, to quantum domain. This is done 
through the extension of the phase space by introducing the momenta tt^ and iTp conjugate to q and p, respectively 
and the subsequent extensions of the lagrangians, hamiltonians, Poisson's brackets, etc. In its full generality, the 
theory describes the dynamics of the quantum ensembles. Its pure state case is reducible to the conventional quantum 
mechanics in q- or p-spaces, however, with a definite rule for ordering of the factors of non commuting operators. 
The latter feature is a direct consequence of the independence of q and p that is maintained at all stages of the 
formalism. Simple rules for assigning an operator F{q,T:q) in q-space to a function F{q,p) in phase space and vice 
versa are prescribed. Extended canonical transformations enable one to go from one extended phase space to another. 
Correspondingly the associated unitary or similarity transformations in the function space enable one to generate 
further state functions from a given one. This unifying feature of the theory makes the comparison of the various 
functions existing in the literature possible and transparent. 

To demonstrate the simplicity and the power of the formalism certain examples are worked out. Treatment of 
Bloch's equation, partition functions for sirnple harmonic and linear potentials, and the mathematical lemmas of the 
Appendix serve this end. Nasiri & Safari Q| and Razavi have found the presented formalism of considerable 
assistance in their study of dissipative quantum systems. 

Acknowledgments: We thank M. R. H. Khajchpour for fruitful discussions on all aspects of the paper and R. F. 
O'ConncU for valuable comments on our treatment of Bloch's equation. 



• Evaluation o{ Xa ~ UaX ■ 

By three Fourier and inverse Fourier transformations we convert x = "0(9)0* (p) exjp{—ipq/h) into the following 



By this provision we have moved both q and p variables to the exponent. Next we expand Ua of Eq H28(l in 
power series, operate by it on the q and p exponents and arrive at 



APPENDIX A 



forms 




(Al) 



k—co 

E 

k=0 




(A2) 



Using Eqs (Al) and (A2) in the expression Xa = C/qX, and inverting (/)'s back to -^'s gives Eq (|29|l 




(A3) 
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As mentioned earlier, for a = 1/2 one recovers Wigner's standard state functions. Q.E.D. 

• Evolution equation for Xa 

To prove Eq (|30|l . it is sufficient to evaluate H.a — UaHUl, where H is the extended hamiltonian of Eq (|10|l . It is 
easy to show that Q — UaqU^ = q — airp and P = UaP = p — aiTq, which is the essence of the transformations 
of Eq (123. We also note that 

Uaq"p"''Ul ^{q- aiTp^ip - aTr,)™ = {p - aTTgY'^q - airX ■ (A4) 

We leave it to the reader to verify Eq (A4) for him/herself for some small n and m. It is needless to say that 
[Q,P] = 0, because the transformation is unitary. With these provisions one finds 

Tia ^ H [q - aTTp,p+ (1 - a)7rg] - H [q - {1 - a)'Kp,p~ qtt,)] . (A5) 

Expansion of the hamiltonian about gives 

2m oq'^ m oq n! oq" op^ 

The Wigner case is for a = 1/2 

p d v-^ 1 ^ Q^ri+ly ^2n+l 

= 5^ + (2„ + 1)! ^2z^'"^' 9g2„+i q 2„+i • (^7) 

n— ^ ^ 

Q.E.D. 

• Ordering rule in a-representation, proof of Equation 
With Eqs (jHH) and (A3) we have 

p + OC 

<g>'">„=/ q^ilj(q-aT)i:*{q+{l-a)T)p'^e'P''/^'dqdpdT. (A8) 



Writing as (ifi)™^-^, integrating by parts m times with respect to r frees the integrand from the p^ factor. 
Then integration with respect to p gives 5{t). Thus 

9"(-*^)"7n;r [^(9 - + (1 - ")^)] '5(r)dqdr. (A9) 

-oo 

Next we substitute ^ by ^ with appropriate adjustments and carry out integrations by parts over q where 
ever necessary to free ip* and arrive at 



f + OO 



i,{q)dq. (AlO) 



' r=0 \ r 

The expression in the integrand is the desired ordering of Eq 1)33(1 , corresponding to q^p"^ in a-representation. 
For a = 1/2 one recovers Weyl's ordering 



1 ™ /m\ 

(2)" El ^ 1 



{t^T > J TT/g^TT,"-'-. (All) 
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The combination of a = and 1, corresponding to averaging by % + is the symmetric ordering of Eq 
(Eg.Q.E.D. 



• The product rule, proof of Eqs (|38|l and (|39() 

The phase space fmiction corresponding to the product of two operators F = AB, by the definition of Eq 136|l . 



IS 



F{q,p) =< q\AB\p > e^'P'^^" = J < q\A\p' X p'\q' X q'\B\p > e~'P'i'^dqdp 

A{q,p')B{q',p)exp ~ 'f^'^^^ dqdp, (A12) 

where by Eq (|35|l . we have substituted < glAjp' >= A{q,p') exp{ip'q/h) and similarly for < q'\B\p >. With 
further change of variables q' — q = q" and p' — p = p" , we obtain 

Fiq,p) = / Aiqy+p)Biq" + q,p)e-^"P" /^dq"dp" = £ ^''f f ^^ ff\ (A13) 

where we have Taylor-expanded A(q,p+p") and B{q + q" ,p) about {q,p) and carried out the required integration 
by parts. 

To deduce Eq H39|) . we first Fourier-transform A{q,p) to a{p', q') and B{q,p) to b{p" , q") in Eq (A13) and carry 
out the necessary differentiations: 

F{q.p) = 2^ -^Q-^ J <^{P,q) exp dp dq—Jb{p,q) exp dq dp 



, , -ip'q + iq'p „ „ -ip"q + iq"p -iq'p" ih, i ^ i , n , ^ata^ 
i[p ,q ) exp b[p , q ) exp e ^ ' dq dp dq dp . 1^.14) 

Next we operate on Eq (A14) by a Taylor-expanded form of Ua as in Eq (A2) and perform the required 
differentiations: 

F^iq,p) = U^Fiq,p) ^ J e-^'^V^exp '^^^^ + '^'^ ' "^^^^^ + aiq\p') 

xe-^'V'/'i exp + W wy,^ q")dq'dp'dq"dp". (A15) 

The exponentials preceding a(p', g') can be written as 

^had^/dpdq -V(g + ag") + iq'jp - (1 - aK) 
exp ^ , 

where the first factor is simply Ua{q,p) independent of the integration variables {q', p', q", p"). With this 
provision integrations over q' and p' can now be carried out and a(p', q') inverse-Fourier transformed. One finds 

Fa{q,p) = / {UcA[q + aq",p- (1 - a)p"]} e'"^"^"/^'e '""'i-^""'' q")dq" dp" . (A16) 
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We again apply the same trick. To the left of the right most exponential we replace, every where, g" by 
{—ihd/dp) and p" by {ihd/dq), perform the inverse Fourier transform of b{p", q") and find 



Fo.{q,p) = Uc,F{q,p) = Ac, 



q + iha—,p— ihail — a) — 
op oq 



Baiq,p) 



E 



{ihY' 



71 = 



d d , . d d 
a— — — (1 — a)- 



Aa{q,p)Ba{q,p). 



dqA dpB dp A dqs 

where d/dpA indicates a differentiation on A{q,p) only, similarly the other differential operators. Q.E.D 



(A17) 
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